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Abstract 



We apply the finite-temperature renormalization-group (RG) to a model 
based on an effective action with a short-range repulsive interaction and a ro- 
tation invariant Fermi surface. The basic quantities of Fermi liquid theory, the 
Landau function and the scattering vertex, are calculated as fixed points of 
the RG flow in terms of the effective action's interaction function. The classic 
derivations of Fermi liquid theory, which apply the Bethe-Salpeter equation 
and amount to summing direct particle- hole ladder diagrams, neglect the zero- 
angle singularity in the exchange particle-hole loop. As a consequence, the 
antisymmetry of the forward scattering vertex is not guaranteed and the am- 
plitude sum rule must be imposed by hand on the components of the Landau 
function. We show that the strong interference of the direct and exchange 
processes of particle-hole scattering near zero angle invalidates the ladder ap- 
proximation in this region, resulting in temperature-dependent narrow-angle 
anomalies in the Landau function and scattering vertex. In this RG approach 
the Pauli principle is automatically satisfied. The consequences of the RG 
corrections on Fermi liquid theory are discussed. In particular, we show that 
the amplitude sum rule is not valid. 
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I. INTRODUCTION 



In 1956-1957 L.D. Landau formulated his theory of Fermi hquids.S The original phe- 
nomenological formulation of this theory is based on an expansion near the ground state of 
the energy functional in terms of variations of the distribution function (bosonic variablesL 
Later, Pomeranchuk derived the thermodynamic stability copditions for this functionalB 
Much effort has been dedicated, including by Landau himself,El to vindicate some intuitive 
assumptions of Landau and elucidate the foundations of the phenomenological Fermi Liq- 
uid Theory (FLT). The field-theoretic interpretation of the Landau FLT has reformulated 
the key notions and basic results ofrtl^ phenomenological theory entirely in terms of the 
fermionic Green functions technique.ou The demonstration of the equivalence of the field- 
theoretic results obtained from the solution of the Bethe-Salpeter equation with the results 
obtained from the functional expansion and from the Boltzmann transport equation de- 
scribing the collective modes, has become a textbook topic.&1a The field-theoretic approach 
provided not only a solid basis to phenomenology, but also a potentially efficient method to 
calculate the phenomenological parameters of FLT from first principles. 

Current interest in non-Fermi Liquids in d > 1 inspired a new wave of efforts aimed at 
clarifying the foundations of the Landau FLT and the mechanisms of its breakdown. Let us 
mention only two approaches, which can be seen as sophisticated modern counterparts of 
the two classic formulations of the Landau FLT. A bosonized treatment of Fermi liquids has 
recently been-developed" in the framework of Haldane's formulation of higher-dimensional 
bosonization.lHl At about the same time, the Renormalization Group (RG) technique has 
been applied to interacting fermions in c? > 1 with models based on fermionic field effective 
actions (see Refs |I2]-|2D| and references therein). It both approaches it has been established, 
for models with reasonable fermion-fermion effective interactions, that the Fermi liquid phase 
is stable, whereas adding gauge-field interactions may drive the system towards a Non- Fermi- 
Liquid regime, or may result in a Marginal Fermi Liquid phase, like for composite fermions 
at the half-filled Landau level. I-, 

The RG analysis of FLT presented here and in our previous work,cJ like other such 
analyses already published, starts from a low-energy effective action with a marginal (in 
the RG sense) short-range interaction. However, contrary to other works on the subject, 
our finite-temperature RG approach revealed that, in the Landau channel of nearly forward 
scattering quasiparticles, the effective interaction flows with successive mode eliminations 
towards the Fermi surface, even in the absence of singular or gauge interactions. In other 
words, the action's interaction (coupling function) does not stay as a purely marginal under 
the RG transformation, since its /5-function is not identically zero. From the RG flow 
equations the standard FLT results have been recovered.ll^ 

It was also pointed out, and elaborated later in more detail by one of us together with 
N. Dupuis in Ref. ^0, that the bare interaction function of the low-energy fermion effective 
action cannot be identifled with the Landau interaction function. The latter, along with 
other observable parameters of a Fermi liquid, should be calculated as a flxed point of the 
RG equations.^ Let us briefly give two arguments for this. First, identifying the Landau 
function with the effective action's bare interaction is inconsistent with other standard FLT 
results, due to the role of Fermi statistics. Indeed, in a stable Fermi liquid, the well- 
known relationship between components of the scattering amplitude (F;) and of the Landau 
interaction function (-F)), i.e., F; = Fi/{1 + Fi), cannot satisfy the Pauli principle for the 
amplitude (the amplitude sum rule) if F has the symmetry properties of the action's bare 
interaction. (For the explanation of this point see Sect. [V| below). Second, identifying the 
Landau function with the bare interaction is inconsistent with the low-energy effective action 
method itself, in the way it is applied to condensed matter problems. Namely, at the starting 
point of the analysis, the bare parameters-of the effective action, including the interaction, 
are regular functions of their variables. It is known, however, that this is not the case 
even for parameters of a normal Fermi liquid. For instance, the scattering amplitude and 
the Landau function are two distinct limits of the four-point vertex in the Landau channel 
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when energy-momentum transfer goes to zero. The non-analyticity of the forward scattering 
vertex appears in its dependence both on the small energy-momentum transfer and, due to 
the antisymmetry (crossing symmetry), on the small angles between incoming (outgoing) 
particles lying near the Fermi surface. This contradiction becomes fiajgrant if one couples the 
fermionic action with gauge fields since, as shown by other methods,l=il the Landau function 
for the marginal Fermi Liquid of composite fermions at the half-filled Landau level develops 
a delta-function singularity in the forward direction (6 = 0). Such behavior of the Landau 
function is related to the divergence of the quasiparticle's effectivg,mass, according to the 
theory of Halperin, Lee and Read for the half- filled Landau leveled (see also Ref |2^) . So, 
coming back to our arguments, the Landau function cannot be a regular interaction in the 
effective action at the starting point of the RG analysis. 

The aim of the present study is twofold. Once the classic FLT results have been recovered 
by the RG approach,t3 the latter would loose its appeal if it did not provide a constructive 
method for calculating the Fermi liquid's parameters. This is especially important goal 
in the long-term prospective of applying this powerful method to more complex strongly 
correlated fermion systems. In this work we explicitly derive the Landau function and the 
forward scattering vertex from the short-range effective bare interaction. We do it in the 
one-loop RG approximation which takes into account contributions of the direct [ZS) and 
exchange {ZS') graphs. This enables us to reveal singular features of the Landau function 
and scattering vertex in the forward direction (^ = 0). 

An equally important goal of this work is to resolve the old problem of FLT with the 
Pauli principle. In its treatment of FLT, the field-theoretic approach encountered a very 
subtle problem caused by Fermi statistics of one-particle excitations and by the necessity 
to provide both stability |£Qp,the Fermi liquid and a solution for the two-particle vertex that 
meets the Pauli principle.cillZl The problem was "settled" by imposing the amplitude sum rule 
on the components of the Landau quasiparticle's interaction function. The phenomenological 
FLT is spared from this problem partially by the way it is formulated, partially because it 
says nothing about the quasiparticle scattering amplitudes. (A detailed discussion of this 
problem, which lies at the heart of the present study, is postponed until SecjV], where it will 
be put in contact with the present RG approach.) The same problem arose m our previous 
workE3 in the form of a "naturalness problem"Ej of the effective action: the effective action 
had to be "fine tuned" in order for the scattering amplitude to meet the Pauli principle. 
We will show that if quantum interference of the direct and exchange processes is taken into 
account, this problem is eliminated in a natural manner. 

The paper is organized as follows. Sections |I| and |ITT| are introductory: we define the 
effective action of the model and the coupling functions (the bare interaction) and vertices to 
be calculated in the Landau interaction channel. In Section which is rather technical, the 
one-loop RG equations for the two-dimensional case are derived. Section explains some of 
the weak points of the standard FLT results and argues for their partial revision. In Section 
[VI| we give a numerical and a ppro ximate analytical solution of the coupled RG equations for 
spinless fermions. In Section |V11| we present and discuss our results for the Landau function 
and the scattering vertex calculated at different temperatures. In Section |V111| we relate 
this study to the standard treatment of Fermi liquid Theory. The consequences of the RG 
corrections on FLT results are discussed. 



II. THE MODEL 



We apply the Wilson-Kadanoff renormalization scheme in the framework developed ear- 
lier for a model with S'f/(A^)-invariant short-range effective interaction and rotation invariant 
Fermi surface in spatial dimensions < 3 at finite temperature In order to make the dis- 
cussion as clear as possible, we concentrate in this work on 2D spinless (A^ = 1) fermions. 
This simple model has nevertheless all the necessary qualities to illustrate our key points 
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and to demonstrate the new features brought by the RG analysis of a Fermi hquid. In this 
case the RG equations take their simplest form, since only the antisymmetric momentum- 
frequency dependent parts of the interaction and vertices are present (they were labeled by 
A in Ref. ^ . 

The partition function in terms of Grassmann variables is given by the path integral 



Z = j VipVifj e^o+^-' (1) 
wherein the free part of the effective actionliHii is 

So= [ ^(l)K+/i-e(Ki)]^(l) . (2) 
-'(1) 



'(1) 

We introduced the following notation: 



(i) 



f3j (27r) 



(i) = (K,,^,), (3b) 

where (3 is the inverse temperature, /i the chemical potential, cjj the fermion Matsubara 
frequencies. We set ks = h = 1. The interacting part of the action is 

Sint = -] [ V'(l)V'(2)^(3)V^(4)r^«(l,2;3,4)/5(27r)25(l + 2-3-4) , (4) 

^ ^(1,2,3,4) 

where S{- ■ ■) stands for a Dirac delta function for the momenta and a Kronecker delta for 
the Matsubara frequencies. The function is antisymmetric under the exchange (1^2) 
and (3 4). The bare cutoff Aq of the action is introduced such that each vector Kj in 
the effective action lies in a shell of thickness 2Ao around the Fermi surface. We denote 
this shell, i.e., the support of the effective action in the momentum space, as C\^^. The 
Matsubara frequencies are allowed to run over all available values. We presume that the 
density of particles in the system is kept fixed. 

The one-particle excitations are linearized near the Fermi surface, and therefore the bare 
one-particle Green's function for the free part of action 5*0 is: 

Gq "^(Ki, tui) = itoi + /i — e(Ki) ^ iuji — vf{Ki — Kp) = iuji — vpki . (5) 
In the integration measure only the relevant part is kept: 

dK = {Kf + k)dkde ^ Kp / dkdO (6) 

J-Ao Jo J-Ao Jo 

The temperature T is restricted by the condition 

T < vpAo . (7) 

The relevant physical information can be obtained by studying interactions of particles 
scattering with small momentum and energy transfer (we call it the Landau channel), and 
those with nearly opposite incoming (outgoing) momenta (the BCS channel). Since we are 
interested in the repulsive case, we presume that stability conditions against Cooper pairing 
are fulfilled, and we concentrate on the Landau channel. 
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III. COUPLING FUNCTIONS AND VERTICES IN THE LANDAU CHANNEL 



Let us clarify the meaning of the quantities entering the effective action. Consider the 
vertex function r(l, 2; 3,4), constructed from the connected two-particle Green's function 
6*2(1) 2; 3, 4) = — {ip{l)ilj{2)ilj{3)il!{4)) c by amputation of the external legs. Here (...) means 
an average with the effective action (0,^) which contains only "slow" modes, lying in the 
support C\^. Once auxiliary source fields (with momenta inside the shell Cj^^) coupled 

to the action's Grassmann fields {ip, ip} have been introduced, such connected n-particle 
Green's functions can be defined as functional derivatives of the source-dependent generating 
functional.@ At tree-level, 1(1, 2; 3, 4) |tree = r^°(l, 2; 3, 4). The bare vertex 1^° (in the 
sense of the effective action (0,§)) can be defined in the same fashion as F, with the difference 
that r^" is the result of averaging over the "fast" modes (those outside C\^) with the 

microscopic action. Contrary to F, the vertex F"^" is not a physical observable, since it is 
not the result of an integration over all degrees of freedom. 

Taking into account momentum and frequency conservation, we use the following nota- 
tion for the nearly forward scattering vertex: 

F(1,2;1 + Q,2-Q) = F(1,2;Q) , (8) 

with the transfer vector 

Q = 3-1 = {Q,Q) (9) 

such that Q ^ Kp is a bosonic Matsubara frequency). We write the momentum K, as 
Kj = + kj where lies on the Fermi surface and kj (|ki| < Aq) is normal to the Fermi 
surface at the point K^. 

In order to calculate physical quantities, we must perform an average with the effective 
action (0,^), i.e., we must integrate out the "slow" modes, which lie inside C^^, in the corre- 
sponding path integrals. This is done in Wilson's RG approach by successively integrating 
the high-energy modes in C^^, i.e., by progressively reducing the momentum cutoff from Aq 
to zero. We define a RG flow parameter t such that the cutoff at an intermediate step is 
A{t) = AqC"*. Integrating over the modes located between the cutoffs A{t) and A(t + dt), 
a recursion relation (in the form of a differential equation) can be found for the various 
parameters of the action. This equation (or set of equations) is then solved from t = to 
t — *■ cxD and this yields the fixed-point value of the parameters of the action. The physical 
quantities are then obtained from these parameters, e.g. by functional differentiation if they 
are source fields. 

A considerable simplification of this problem comes from the scaling analysis of the low- 
energy effective action using the smallness of the scale A/Kp-^ A tree-level analysis shows 
that the only part of the coupling function F^ which is not irrelevant couples two incoming 
and two outgoing particles with the same pairs of momenta (K]?, K|,) lying on the Fermi 
surface. The dependence of the coupling function on kj and on the frequencies Ui is irrelevant 
and can be omitted. When the initial cutoff Aq satisfies condition (0), we can unambiguously 
define a bare coupling function which depends only on the angle between the incoming (or 
outgoing) momenta. This bare coupling function is given by the vertex F'^o(l, 2; Q) in the 
zero transfer limit (Q = 0) where the two external momenta are put on the Fermi surface 
and the external frequencies are c^min = irT (the latter will be dropped from now on). 

f/(K^, K|) ^ ^upT^^iK'p, Kl; 0) , (10) 

where = Kf/t^vf is the free density of states at the Fermi level. Each vector may 
be specified by a plane polar angle 6'j. The function U is an even function of the relative 
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angle 9i2 between and_K|,. The only remnant of the antisymmetry of T^" (the Pauli 
principle) is the condition:ll2l 

U{0) = . (11) 

As shown earlier,^! the tree-level picture becomes more complicated when we carry out 
the mode elimination inside C\. It turns out that simply discarding the frequency depen- 
dence of and identifying the momenta K.p ^ K|,, K|, ^ is an ill-defined procedure 
when the running cutoff A becomes of the order of the temperature [vpA ~ T). The ambi- 
guity arises when calculating the one loop-contribution from, say, the-.ZvS, graph, since this 
contribution is not an analytic function of the transfer Q at Q = oBlro To describe cor- 
rectly the parameters of the Fermi liquid, one should retain the dependence of the coupling 
function T^{K.]p, K|.; Q) on the energy-momentum transfer Q. Retaining this Q-dependence 
allows the calculation of response functions or collective modes of the Fermi liquid.£^ For 
the purpose of the present study we define two coupling functions (F*^ and F^), depending 
on the order in which the limits of zero momentum- (Q) and energy-transfer (Cl) are taken: 

F«(^i2) = lim [f(^^i2, Q) 1 , (12a) 



n=o. 



'12 J 



Q=oJ 



) = lim T{9i2, Q) (12b) 



We use dimensionless vertices, by including in their definition the factor ^i'f, like in Eq. (PH]). 

The functions F'^'^(^) are even functions of the angle 6. We will not explicitly indicate their 
dependence on the cutoff A, unless necessary. We will indiscriminately call these functions 
(running) vertices. 

Let us summarize: The effective action is defined on the support with the bare 

coupling function F^°, which is presumably an analytic function of its variables and is 
marginal at tree level. While performing the mode elimination within C^^, we need to 

calculate the flow of the two vertices F'^ and F^. The bare coupling F'^" has an unambiguous 
meaning only as the common initial point of the RG flow trajectories of F*^ and F^. The 
fixed point values F*^* = T'^{t = oo) and F^* = T^(t = oo) are physical observables: the 
first one is the Q-limit of the vertex F (as defined at the beginning of this section) and is 
the scattering amplitude of quasiparticles with all four external momenta lying on the Fermi 
surface. The second one i^ the unphysical limit (f2-limit) of the vertex F and is identified 
with the Landau function.^ 



IV. THE RG EQUATIONS IN THE LANDAU CHANNEL 

There are three Feynman diagrams contributing to the RG flow at the one-loop level (see 
Fig. 1), denoted ZS (zero sound), ZS' (Peierls), and BCS. The BCS graph contribution 
preserves the antisymmetry of the vertex, while those of the ZS and ZS' graphs separately 
do not: only their combined contribution [ZS + ZS') is antisymmetric under exchange of 
incoming (or outgoing) particles. To respect the Pauli principle, it is therefore necessary to 
take into account both the ZS and ZS' contributions to the RG flow. In this work we discard 
the symmetry-preserving contribution of the BCS graph to the RG flow of the vertices in 
the Landau channel. Thus, we leave out the interference near 9 of the Landau channel 
with the BCS channel, which leads to the Kohn-Luttinger effect .li^ 
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The formal analytic expression of the ZS graph is 



ZS 



r(l, 5; 1 + Q, 5 - Q)r(5 - Q, 2; 5, 2 - Q)G(5)G(5 - Q) 



(13) 



(5) 



wherein the transfer vector Q is given by (^. To calculate the contribution of this graph to 
the RG flow of F*^ and F^, we only need to keep the dependence on the momenta and on 
the transfer Q in the vertices on the r.h.s. of (|l^). Momentum and energy conservation is 
already taken into account in ([T3|). The phase space restrictions are satisfied automatically 
for any K5 G in the limit Q ^ 0. When Ki and K2 lie on the Fermi surface and Q — 0, 
the r.h.s. of (O) contains both vertices of type ([T^ ) with Kp running freely around the 
Fermi surface during the angular integration. Thus, for this graph, all the phase space is 
available for integration. The summation over of the Green's functions product on the 
r.h.s. of ([I3|) when Q — > gives zero in the f2-limit, and thus 



dt 







(14) 



ZS 



The Q-limit of the same product gives a factor \f3 cosh, '^{jdvph^/'^), and accordinglyS 



dt 



4 



ZS cosh {I3r) 



(15) 



where we introduced a dimensionless temperature flow parameter: 



Pnit) = -VF(3A{t) ■ 
We now turn our attention to the ZS' graph. Its analytic form is 

ZS' = [ F(l,5;l + Q',5- Q')r(5- Q',2;5,2- Q')G'(5)G(5- Q') 



(16) 



(17) 



wherein Q' = 2 — 1 — Q can be thought of as an "effective" transfer vector for this graph. 
For IK2 — Ki| 7^ the limit Q — of the r.h.s. of ( p!7D is single- valued and equivalent to 
the Q-limit. BJ The Green's functions contribution to this graph is 



Wg{5)G{5-Q'] 



tanh 



/i) 



tanh 



(e(K. - K. 



Ki)-/i) 



Q=0 



Ki; 



(18) 



If l^^i — ^^2! ^ T/vpKp the r.h.s. of Eg . (p!8|) becomes —j/3 cosh '^{(3vFk^/2). The calculation 
of the ZS' contribution is more subtle, since even in the zero-transfer limit Q — (in any 
order), the vector Q'|q^o = K2 — Ki is free to take any modulus in the interval [0, 2Kp\ as 
the angle Oi — 62 varies. A large Q' kicks the vertex momenta on the r.h.s. of ([17| ) outside 
of C|, even if K5 G C|. In such cases the contribution of the ZS' graph is cut off, except 
for special positions of the vector K5 running over the Fermi surface. Thus, for an arbitrary 
angle 61 —62, not all the phase space is available for integration. 

To understand where this elimination of the ZS' contribution comes from, we must keep 
in mind that our effective action has support C\ in momentum space. Let us consider the 
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ZS' graph (see Fig. 1) when all external momenta satisfy momentum conservation and lie 
in C\. It suffices then to check whether the internal momenta (K5 and K5 — Q') lie in C\ 
when K|n runs around the Fermi surface during the integration. From Fig. 2A we see that 
if IQ'I > 2A, the loop momenta lie both in C\ only at special values of K|. (the shaded 
regions), i.e., only small fragments of phase space are available for integration. At smaller 
Q' (cf. Fig. 2B) these intersections form a connected region and K|. is free to run around the 
Fermi surface. If we completely neglect the ZS' graph when the intersection is disconnected 
(in Fig. 2A), the contribution of this graph to the RG flow at |Q'| < 2A is calculated in the 
same way as that of the ZS graph. Since |Ki| = IK2I = Kp and Q'Iq^o = K2 — Ki, the 
condition |Q'| < 2A is equivalent to the condition | sin((6'i — 6*2)72)1 < kjKp for the angle 
between Ki and K2. 

Taking into account both the contributions of the ZS and ZS' graphs, the RG equations 
for r*^'^ can be written in implicit forni-x^ 



dt 



dt 



zs 



'~dr 



(19) 



ZS' 



dT^ 



dT^ 



ZS' 



~dr 



ZS' 



dt dt 

Summing up all formulas, we obtain the following system of RG equations: 

9rQ(20) (3r 



dt 



(20) 



(21a) 



gr"(20) 

di 



— r«(0 - 9)^^(9 + 0)y(0, 9; Pn) • (21b) 



To simplify those formulas we parametrized the angular dependence of the vertices in 
Eqs. d^) by the angle between K], and (K], + K^,), |0| e [0,7r/2]. The small ZS' 
contribution coming from |sin0| > A{t)/Kp (Fig. 2A) was neglected, which is accounted 
for by the Heaviside step function B, wherein 9c = 8iTcsm{A{t) / Kp) . We also defined the 
function 



Y{(P,9;Pr) = 



sinh(2/5Q' 



cosh(2/3ij) + cosh(2/3Q/ 



(22) 



Pqi = Ppsm9 sincj), Pp = f3vpKp . 
which arises in the calculation of the ZS' contribution (O). Notice that 



lim F(0,0;/3h) 



1 



cosh Pr 



(23) 



(24) 



From Eqs. (^T]a,^^) we see that at small angles (|0| < T/vpKp) there is a strong inter- 
ference between the ZS and ZS' contributions. This interference depletes the RG flow of 
r'5(0) at small angles. Moreover, at = the flow is exactly zero, for the two contributions 
have the same thermal factor cosh~^(/5ij):il 
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(25) 



The initial conditions for the flow equations (^Tj) are: 

r«(0,t = o) = r^(0,t = o) = f/(0) . (26) 

Recall that the fixed points F*^* and F^* of the vertices F*^ and F^ are the forward scattering 
vertex and the Landau interaction function, respectively. From Eqs. (|25| , p6| , p!T|) we conclude 
that the RG equations for the forward scattering vertex preserve the raun principle at any 
point of the RG flow trajectory 

r'3(0 = 0,t) = , Vt, (27) 

while the "uncompensated" RG flow generated by the ZS' graph drives the vertex F^ to 
a fixed point value (the Landau function), which does not satisfy the Pauli principle, i.e., 
T^*{(f) = 0) ^ 0. 



V. DEFICIENCIES OF THE DECOUPLED APPROXIMATIONS IN THE 

LANDAU CHANNEL 

Before finding a solution (exact or approximate) to the flow equations (PTT) which fully 



takes into account the coupling of F*^ and F^, we will comment on approximate solutions in 
which this coupling is neglected. The Landau channel, as defined in this paper, includes, at 
one-loop RG, both the direct {ZS) and exchange {ZS') quasiparticle-quasihole loops with 
a small transfer Q. We will call decoupled any treatment of the Landau channel which does 
not explicitly take into account both the direct and exchange contributions. It is shown 
below that solutions for the forward scattering vertex provided by decoupled methods fail 
to meet the requirements of the Fermi statistics. Tackling the Pauli principle by imposing 
additional constraints on the solutions (sum rules) leads to conceptual difficulties discussed 
below. 

To shorten notation we drop upper labels (Q, f2), and define F {F) as the running vertex 
whose fixed point is the forward scattering vertex (resp. the Landau function). 

Let us first solve the RG equations in the decoupled approximation. If we neglect com- 
pletely the ZS' contribution in_Eqs. ([21| ) and perform a Fourier transformation, we recover 
a familiar system of equations,li2l with its RPA-like solution in which all harmonics are de- 
coupled: 

Bt ' ' ^ ' l + (ro-r)FKro) ^ ' 



^ = =^ F^^\t) = est. (28b) 

We introduced the auxiliary parameter r = tanh/?/? (r G [0,ro]), with tq = tanh/3o and 
Pq = ^vfPAq. Since the temperature in the effective action is restricted by the condition 
(0), we can set tq = 1 for all practical purposes. 

With the initial conditions F;(ro) = Fi{to) = Ui (cf Eq. (P^)), the fixed points of Eqs.( 

are 



9 



(a) rr = ^ (b) F; = Ui, (29) 

with the following stability conditions for the fixed point: 

Ui>-1, V / , (30) 

which are the Stoner criteria well-known from the RPA approach. The bare interaction 
satisfies the Pauli principle (cf. Eq. ([TT|) ) 



J2Ui = 0. (31) 



(=— oo 

If the vertex F is to satisfy the Pauli principle, the condition 



l=~oo 

must be imposed on the r.h.s. of (|29|a). However, it has been known for a long time th 
conditions (152) and 



Indeed, subtracting (11) from (|52D, we find 



are incompatible, unless the stability conditions (|5D|) are broken.c3 

, (33) 



l=—oo 

which cannot be satisfied without violation of (pOD. 

This proves that the antisymmetric bare interaction U cannot be at the same time a 
fixed point of the RG flow and the Landau function, unless the classic FLT formulas are 
unapplicable. The accepted cure to this paradox is to give-|Up the Pauli principle on the 
Landau function, because of the neglected ZS' contribution.c3 In the RG approach, this may 
be accomplished (in the decoupled approximation) by letting the ZS' contribution drive 
the bare interaction U towards the Landau function F* during an earlier stage of mode 
elimination, and then by solving the RG equations (|28|) with F* as a new renormalized 
"bare" interaction.ll3E3 This leads to the well-known relationship between the scattering 
vertex and the Landau function 

rr = ^ (34) 

Because of the ZS' contribution, the Pauli principle does not apply to F* (a), while it is 
enforced on the vertex F* through a sum rule (b): 



(a) J]f;^o (b) ErTF^ = °- 



Ft 

^ T 

1= 

In doing so, the stability conditions (|30|) are modified as follows 

F;>-1, V/, (36) 

i.e., they become Pomeranchuk's stability conditions for the Fermi liquid, originally obtained 
on thermodynamic grounds.^ Such a decoupled RG treatment of the direct and exchange 



loops makes Eqs. (|35D compatible with the conditions (|36|) . 
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However, the sum rule (|35|b) is "unnatural" , in the following sense. The bare interaction 
can in principle be traced rFom a microscopic Hamiltonian. For instance, let us consider 
the spinless extended Hubbard Hamiltonian on a square lattice (with lattice spacing a) 
at low filling, with nearest-neighbor repulsive interaction ([/™). Fourier-transforming and 
antisymmetrizing the interaction, we end up with the following coupling function of the 
microscopic Hamiltonian: U^^^(Ki, K2, K3, K4) ~ - jO^W^ ■ (Ki - K2) ■ (K3 - K4).0 Let 
us choose this interaction as a trial bare dimensionless coupling function: 

U{e,-e,)=Usin'(^^^^) , (37) 

wherein all parameters are hidden within a single coefficient U. The only nonzero Fourier 
components Ui of the interaction are: 

Uo = lu , t/±i = -\u . (38) 



The interaction (|37|) satisfies the Pauli principle 



11,31). The RPA sum rule (p3) imposes 



an additional constraint, which the interaction (p7)~does not satisfy. If we suppose that the 
"improved" results ( p^ , ^5| b) are always true, then, starting from any kind of microscopic 
interaction (e.g., the bare interaction (|37D) and integrating "fast modes" outside the imme- 
diate vicinity of the Fermi surface, we have to end up with a "fine tuned" interaction, for 
any interaction has to be "fine tuned" in order to satisfy (|35|b). The integral of the fiow 
(p7|) (or, equivalently the sum rule (^81) below) is not a fine tuning, since firstly, the bare 
interaction at the initial point can be always antisymmetrized, and, secondly, we have an 
exact cancellation of the RG fiow for the vertex F at zero angle due to direct and exchange 
contributions, thus preserving (p?]). On the contrary, there is no reason for any bare inter- 
action to satisfy (^) at the begmning, nor is there a mechanism to provide the fine tuning 
(P5|b) on other parts of the RG trajectory. 

These difficulties are not specific to the decoupled RG approximatimi, since the latter is 
strictly equivalent to the diagrammatic microscopic derivation of FLTHoQ leading to the same 



results ( p4| , p5| , |36|) . The decoupled RG treatment is equivalent to applying the Bethe-Salpeter 
equation with the particle- hole ZS loop singled out, F being the vertex irreducible in this 
loop. There are no a priori reasons in that approach to demand this vertex to satisfy the 
Pauli principle. The rearrangment of diagram summations in the Bethe-Salpeter equation 
leading to (Q) is based on the assumption that the vertex irreducible in the direct particle- 
hole loop (ZS) is a regular function of its variables, neglecting the zero-angle singularity (at 
T = (£3) in the ZS' loop. As a consequence, the Pauli principle for the scattering vertex F* 
is not guaranteed in the final result and "the amplitude sum rule" (p^) must be imposed 
by hand. The solution (B^ of the Bethe-Salpeter equation is tantamount to the summation 
of the ladder diagrams Duilt up from the ZS loops, wherein the Landau function stands 
as the bare interaction. For this reason, the solution ( p4D we will call the ^HheZ S -ladder 
approximation" in the following. We refer the reader to a paper of A. HewsonH wherein a 
"generalized" Bethe-Salpeter equation for Fermi liquids, which explicitly takes into account 



both the ZS and ZS' loops, is derived. For further discussion on this issue, see also Ref. |2| 



VI. SOLUTION OF THE COUPLED RG EQUATIONS 



A. Exact Numerical Solution 



The coupled integro-differential fiow equations (PTI) may be solved numerically. The 
functions T{6) and F{6) are then defined on a discrete grid of angles, and simple linear 
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interpolation is used to represent them between the grid points. The grid spacing is not 
uniform: it has to be very small near ^ = 0, where the flow is singular, but may be larger 
elsewhere. The RG equations then reduce to a large number of coupled nonlinear differential 
equations, which are solved by a fourth-order Runge-Kutta method with adaptive step-size. 
Typically, a grid of a few hundred points is sufficient (we take advantage of the symmetry of 
the functions). Of course, the numerical solution was checked to be indistinguishable from 
the (exact) RPA solution when the ZS' contribution is discarded. 

An example of solution for the spinless case with the interaction function (E7t) is shown 
on Fig. 3(A), at various temperatures. The interaction function U{9) and the RPA solution 
-pRPAj^^-j are also shown. This solution will be discussed in Sect. |VII| . 



B. Approximate Analytical Solution 

The flow equations (^) may also be solved analytically, albeit only approximately. In 
this section we give the approximate solution for the fixed points F* and F* both in terms 
of Fourier components and in terms of angular variables. (See Eqs. (4^,|^) below.) 

The Fourier transform of Eqs. (^) is 

-F^ + — (39a) 



di^R cosh^/3« m 



R 



d(5. 



R 



n—m,2l—2m 



{PR)riT. 



l-2m 



(39b) 



yn',m'{l3R) = \ r'\(t) / rf^ cos(20n')cos(^m') e(^e-|0|)m^;/?R) • (39c) 







On the plane {4>,0), the function Y{(j),6] /3r) has a maximum on the line 6 = 7r/2, which 
moves from the position (7r/2,7r/2) at the beginning of renormalization procedure (when 
Pr ~ Pf) towards the position (0,7r/2) when approaching the fixed point (/?/? 0). Else- 
where, Y{(j), 6; (3pi) is either quite fiat, or its contribution is eliminated by the cutoff fac- 
tor Q{9c — 101) during the renormalization flow. Therefore, we approximated the function 
Y{(j),6] Pfi) on the plane {(f), 6) by its value on the line (0, 7r/2). This approximation, sim- 
plifying considerably our equations, allows an analytical treatment and a qualitative insight 
harder to find in purely numerical results. The approximate analytical solution of the RG 
equations given below justifies that simplification a posteriori, when compared with the 
direct numerical solution of Eqs. (|2TD. 
The approximate RG equations are: 

5r„ 



E [—r2^S^m-Y^_UPR)]rl (40a) 



dpR ^ Lcosh^/3i? 



d^R 
wherein 



^n-rriPR)-^1 , (40b) 
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2 /•arcsin(2/3H//3F) 

yn{(3R) = - d(t> -; (3r) cos(2n0) . 



(41) 



The key difference between Eqs. (pSD and (TO) is that the former do not generate new har- 
monics since all harmonics are decoupled, whereas the latter couple all harmonics (because 
of the ZS' contribution) in such a way that an infinite number of new harmonics are gen- 
erated by the RG flow, even if only a finite number of harmonics are nonzero at the start. 
For instance, the trial interaction (p3) has only three nonzero components, but according 
to Eqs. ( ^0] ) the fixed points T* and F* will possess an infinite number of them. The gener- 
ation of new harmonics is not an artefact of the approximation which was used to go from 
Eqs. (^) to Eqs. (|40|) , but is a generic consequence of the interference in the Landau channel 
(cf. Eqs. (H). 

Let us start the analysis of Eqs. (|40|) with a heuristic observation. Whereas the compo- 
nent Yq{(3r) is a nonnegative function of [3r, the others (^^(/^r), n > 1) are increasingly 
oscillating functions of (3r when n increases. These oscillations along the whole RG trajec- 
tory [0, /?o] "will effectively decrease the contributions from the harmonics (m ^ n) to the 
flow of F„. Because of this, we expect the diagonal terms (m = n) of Eqs. (|40|) to be more 
important, and this justifies a perturbative approach, in which the nondiagonal terms are 
ignored at zeroth order. Let 7n(/3ij) be the zeroth order solution: 



OPr 



cosh (3r 



In 



(42) 



The solution is 



7n(A 



R) 



1 + 



tanh/3o - tanh/3ij - /o(/3o) + Iq{I3r) 



Un 



(43) 



with 



I^{(3r) ^ / df3'^ r„(/5^) = -- d4> 
Jo ^Pf Jo 

'3 



In 



cosh(/5j:j + Pf sin < 
cosh(/3/j — Ppsini 



In 



cosh (|/3ir sin 0)-| cos(2ri 
cosh(|/9irsin</) 



smi 



(44) 



The fixed point 7* is 



7« 



1 + 



Wo 



(45) 



The integrals In{Po) can be evaluated analytically, since {Pf,Po) ^ 1 according to condition 
(0). In the following we shall need the first two components only: 



Wo) 



A^l 

KpTT 



In 2 + In ■ 



1 + y/l - AI/{2Kf) 
1 + - Al/K'^F 



+ — 2 arcsm — arcsm -— + 

ttV 2Kf KpJ vfKf 



T (In 2) (In 3) 



TT 



(46a) 
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ln2 + ln 



1 + v/1 - K/i:^Kpy 



1 + v/i - K/Kl 



1-AI/KI-J1-AI/{2Kf 



T (In 2) (In 3) 



(46b) 



(The next term in the temperature dependence, omitted in Eqs. (^Bj), is of the order 
{T/vfKf)^)- Treating the off-diagonal terms (n 7^ m) on the r.h.s. of (^0|a) as pertur- 
bations, we obtain the following approximate solution at first order: 



(47a) 



/3r ^ cosh^ 



(47b) 



It is straightforward to check that the solution (|47|a) satisfies the sum rule (i.e., the Pauli 
principle (25|,p^)): 



(48) 



The solution (|47|) can be converted back in terms of the relative angle 6 G [— vr, tt] with a 
little help from Eq. (H): 



T*{e) = U{e)- r E cos(n%^(/?«) 

n=— 00 



cosh /?R 

ft 



+6(^0 - l^^l) 



71 



i/3p.|sin{e/2)| 



c?/?«F(-,-;/3^) ^ cos(n^)7^(/?«) 



(49a) 



F*{e) = rie)+ ^ cos(n^)r^(/5« 

Jo cosh Br J--'_ 



(49b) 



wherein 60 = 2 arcsin(2/9o//3F)- A comparison of Eqs. (^D|a) and (pi|a) shows that - with the 
aforementioned approximation of the angular dependence of the function Y - the approxi- 
mate solution (|47|a) may be obtained by replacin g th e vertex components r„ on the r.h.s. 
of Eq. (p|a) by the "renormalized" RPA ansatz (^31). It would be a mistake, however, to 
conclude that the ZS' diagram contributes only to the third term on the r.h.s. of Eq. (|49| a) 
since the 7„-s partially include its contribution. It is worth noting that Eqs. (|47|b,^9|b) are 
not approximations in the sense of E qs. (^7| a) or (^a), but they are exact relations for F, 
derived from the basic RG equations (|40|) . 



C. Extension of the Effective Action 



In the numerical and analytical results presented in the following sections the initial cutoff 
Aq of the effective action is extended to Kp, i.e., Po = [3f/2. This point should be clarified. 
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Notice first that the ZS contribution is not sensitive to the bandwidth cutoff Aq - provided 
condition (0) is satisfied - since tanh (3o is unity with exponential accuracy. On the other 
hand, the angular cutoff of the ZS' contribution (cf. Eqs. ( pl|,|39l ,0,|4^)) comes from a cutoff 
imposed on the momentum transfer in this graph (cf. Eg. (|17[))^t is 9^ = a,rcsm{2Pji/ Pp) 
(with 2(3r/ (3f = A/Kp) if Aq < Kp, and 6c = Tr/2 otherwiseT?he specific choice /3o = (3p/2 
(Aq = Kp) means that at the initial point of the RG flow the angle is allowed to take all 
values (i.e., the momentum transfer Q' is not cut off), while the bandwidth is extended to 
the full depth of the Fermi sea. It can be checked that the results are not sensitive to the 
choice of a bigger cutoff Aq > Kp, since then not only is the ZS contribution to the flow 
is exponentially small, but that of ZS' as well, until the cutoff decreases to A ~ Kp (this 
was also confirmed by direct numerical tests). The formulas for the approximate analytic 
solution are derived for Aq < Kp. 

Such an extension of the low-energy cutoff to large values is analogous to what is routinely 
done in ID models (e.g., the Tomonaga-Luttinger modeled). In that context, deviations of 
the real excitation spectrum from linearity and the approximated integration measure are 
expected to affect only the numerical values of the renormalized physical parameters. 

Choosing Aq ~ Kp renders the RG fixed points (observables) sensitive only to the two 
independent physical scales present in the model: T and vpKp = 2Ep, and not to the 
arbitrary scale Aq, which divides fast and slow modes. Lowering the running cutoff until 
it reaches some intermediate scale Ax (such that Ax <^ Kp and vpAx ^ T) provides us 
with Ax-dependent parameters for the action. We regard Ax as the scale of the low-energy 
effective action. However, the observable quantities (the fixed points) do not depend on a 
particular choice of Ax- 



VII. ANALYSIS AND DISCUSSION OF THE RG RESULTS 



We will now discuss the main novelties brought by quantum interference in the Landau 
channel and compare with the results of decoupled approximations. The solutions r*(^) 
and F*{6) at different temperatures and for the interaction (|38| ) are shown on Fig.3 (A: 
direct numerical solution of Eqs. (pi]); B: solution (P9|)). For this interaction the sum in 
the second and third terms on the r.h.s. of Eq. (|3a) is 7o(/5i?,) + '^iHPr) cos 6*. The curves 
were calculated for U = 1 (cf. Eq. (|37D), which isfour times smaller than the critical value 
U^^^ = 4 at which the instability appears in the RPA solution (p8|a) for Tl. Comparison of 
the approximate solutions (E^,^) with the direct numerical solution shows good agreement. 

In Fig. 3 the differences oetween the RG solution and the RPA solution (P5|a) are minor 
at large angles, but they become especially striking at small angles 9, where the interference 
between the ZS and ZS' contributions is very strong. The RG solution gives r*(^ = 0) = 
(the Pauli principle), while r^^-^(^ = 0) = — 1/3 for this interaction strength. The Landau 
interaction function F*{6) differs from the bare interaction U{6), and F*{6 = 0) 7^ 0. If the 
ZS' contribution is neglected (the RPA solution (p8|b)), these two quantities coincide. 

An interesting feature of the RG result is the temperature dependence of the vertices 
r*(^) and F*{6). As T decreases, the "beak" of r*{6) in the region of strong interference 
becomes narrower. The characteristic angular width of this "beak" is \6\ ~ T/vpKp. A 
similar narrowing is noticeable in the temperature dependence of F*{6). One can also see 
from the figures a weakening of the interference effect at lower temperatures, for then the RG 
solutions lie closer to the RPA curves, but the distinctioxis between them do not disappear 
as T — »• 0, and the RG never reproduces the RPA result.E3 

In terms of Fourier components this behavior manifests itself in a linear temperature 
dependence of F* and F*. This linearity is found both in the direct numerical solution of 
Eqs. (|21|), and from the solution of Eqs.( |41|j43 

revealed analytically. Integrating by parts and using Eq. (^), we can rewrite Eq. ([47|a) at 
the fixed point as 



mm. 



This temperature dependence can be 
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r; = 7: + E in-umi - 2 ^ df3n (/^- ^ ^ ^ 



cosh 



YoiPn) iUPr) . (50) 



The leading term on the r.h.s. of Eq. (pOD is 7*. Using then Eq. 
n = 0, 1: 



la), we obtain for 



wherein 



r:(T) ^ 7:(T) ^ 7:(o) + 



7:(o) 



T (In2)(ln3) , 

-bnW] , (^^ = 0,1) 



71 



1 + 



l-/o(/?o 



T=0J 



Un 



(51) 



(52) 



For the interaction (^) t/„ = and so 7* = for n > 1. Thus, the higher harmonics r*>^, 
are entirely generated by the RG flow. To leading order, we obtain from Eq. (|50|) : 



(53) 



This component also has a linear temperature dependence, according to Eqs. (^6[). To 
estimate the components of the Landau function, we first rewrite Eq. (0b) in another, 
equivalent form (cf. Eqs. (|40|)): 



00 

JPr 



(54) 



m=— 00 



Proceeding in the same fashion as above, we obtain the linear temperature-dependent com- 
ponents F*: 



F: ^Un + Wo)Ul + {\n-l\ + l)h{(3o)U\n-iv (^ = 0' 1) 



(55a) 



F;^h{l3^)Ul. (55b) 

We should emphasize that simple formulas like (|5T| , ^j55D serve only to illustrate how the 
temperature dependence comes about, and give omy the order of magnitude of the higher 
harmonics (n > 2). The latter should rather be calculated numerically. The temperature 
dependence of the lowest harmonics (e.g., Fq and F^) does not seem to be a relevant issue 
in the calculation of quantities such as the compressibility, effective mass and heat capacity, 
since, in the total ZS' contribution, the temperature corrections, of the order of T/vfKf, 
are very small in comparison with the main corrections of order Ao/Kp. As a consequence, 
the actual values of the lowest harmonics vary within a few percent at most, even in the 
entire temperature interval < T/vpKp < 0.1 (the maximum temperature studied is really 
high: T = 0.2Ef). 

The temperature dependence is more pertinent as a "collective" effect of the higher 
harmonics generated by the RG flow. Let us explain this point with the example of the 
interaction (PH). The "improved" RPA ansatz (ffil) renormalizes the bare components Un 



into 7„ (n = (J,±l). The latter form almost perfectly the function r*(6'), except at small 
angles. For those three components 7„ the sum rule (EB|) is less violated than for the "pure" 
RPA components ([25|a). The generation of the new harmonics by the second term on the 
r.h.s. of Eq. (P7|a) gives "a final touch" to the curve r*(^), resulting mostly in the formation 
of a temperature-dependent feature near 6* = 0. The actual calculation of the components F* 
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showed that, in order to obtain with acceptable accuracy the right form of r*(^) provided by 
Eq. (|49|a) via the Fourier transformation of Eq. (^7|a), at least A^max ~ vpKp/T components 
are necessary. So, the lower the temperature is, the more harmonics are needed for the 
formation of the vertex T*{9). The same conclusion can be drawn from a numerical solution 
of the equations, but since it is carried out in terms of angles on a discrete grid, a reliable 
calculation of higher harmonics is difficult. 

Another physical consequence of the quantum interference in the Landau channel is the 
increased robustness of the system against instabilities induced by strong interactions. Even 
from the approximate solution (^7[), we see that the maximum interaction strength allowed 
is now larger than the one provided by the RPA solution (cf. (^^D|)). From Eq. we 
obtain the stability conditions for the approximate solution (^71): Ui > —\1 — Jo(/?o)] ; V / 
with < Io{(3q) < 1 according to (^a). Since /o(/?o) grows with temperature, larger values 
of \Ui\ are allowed as T increases: the higher the temperature, the more stable the system is, 
as it should be from physical grounds. At the optimal choice of the initial cutoff (Aq = Kp)-, 
loiPo) grows from 0.255 at T = to 0.27 at T/vpKp = 0.1. This value of temperature is the 
largest we can try without violating the condition of applicability of our model ^). Thus, 
within this approximate solution, the effect of interference increases the critical coupling 
by 40% compared the the RPA critical value (pO]). Since we are retaining only two one- 
loop diagrams, linearized excitation spectrum ana integration measure, we cannot be more 
conclusive on the role of the modes deep into the Fermi sea in screening a microscopic 
interaction of arbitrary strength, and in stabilizing the Fermi liquid phase. 



VIII. CONTACT WITH THE LANDAU FLT AND DISCUSSION 

In this section we|^xplain how the present RG theory is related to the standard results 
of the Landau FLT.Iila This will akQ_.allow us to relate this study to previous work on this 
RG approach to the Fermi liquid. t^'Ey 

It is important to notice that the two contributions to the RG flow, coming from the ZS 
and ZS' graphs, behave quite differently as the flow parameter runs from /^o ^ 1 towards 
(3fi = 0. At large I3r the ZS contribution to the flow, which gives the term proportional to 
cosh~^ on the r.h.s. of Eq. (|5PD, is virtually negligible, up to ~ 1. On this part of the 
RG trajectory, the main contribution to the renormalization of F and F comes from the ZS' 
graph. On the other hand, closer to the fixed point {i3r < 1), the ZS contribution grows 
since cosh~^/5ij ~ 1 for all harmonics, while Yn{PR) decreases for the lower-order harmonics. 
At /3r < 1: 

ynipR) ~ — sm — — . (56) 



vrn P 



F 



Using the approximated form (41) is justified here, since at (3b. ^ 1 there is no difference 



between the exact form of the RG equations (|39|) and Eqs. (^). Indeed, when -C 1, the 



largest allowed (f) is roughly 2[3r/ (3f, so in Eq.(p2|) max|/?Q/| f» 2[3r <^ 1 and the limit (^if) 
of the function Y can be taken. The Kronecker delta appearing after the integration over 
9 removes one summation, and we recover exactly Eqs. (40]) with l^(/5/j) given by (^6|). It 
should be also kept in mind that the ZS' flow is localizedwithin the angle |0| ~ 2(3rJpf- 

Such different behavior of the two contributions {ZS and ZS') to the total RG flow 
explains^aii^ approximations based on the decoupling of these two contributions (RPA, ZS- 
ladderB'LHlHlEl) are reasonable. To clarify to what extent the standard results of FLT ( P^|55[ ) 
can be corroborated by RG, we will make a two-step approximation of our RG equations. 
In doing so we will follow exactly the "recipe" of the ZS'-ladder approximation discussed 
in Sec.V, but now we can check each step by direct comparison with the RG solution of 
Eqs. dill). 
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In the first step we neglect the contribution of the ZS graph above an intermediate flow 
parameter f3x- As one can see from the RG equations (|39|) , this removes the exponentially 
small difference betwefin r„(/3ij) and -Fn(/3fl) at {3^ > /5x~This approximation is asymptoti- 
cally exact as T O.EZI Neglecting, in the second stage of this approximation, the ZS' flow 
for Pji < Px, localized by that time within the angle 9x = 2 a,Tcsm{2Px / Pf) , we recover the 
exactly solvable equations (|28| ) with the new initial point Pr = Px, instead of Pr = Po- Then 
according to Eqs. (0), = Fn{Px) is the (approximate) fixed point value of the Landau 
function, while Vn{PR) flows towards the (approximate) fixed point from the new bare 
value V-^ = Tn{Px) = ■ This second step of approximation violates the Pauli principle, 
no matter how close we are to the Fermi surface (cf. Eq.(^) and Ref. |2^). Afterwards the 
theory says nothing about the values of the functions r(^) and F{9) inside the interval 29 x 
and, of course, there are no more correlations between these functions. 

To preserve the correct zero-temperature limit and to minimize the angle within which 
the approximation gives completely wrong results for F* and F*, the intermediate cutoff Ax 
corresponding to Px = vfAx/2T should be chosen such that tanhPx ~ 1 (cf. Eqs. (p8| , p9[) 
and Ref. [19|) and 2Px/Pf = ^x/Kp -C 1. Summing up what is said above, we obtain: 

pX jpX 

" l + tanh(/5x)r^ l + F^ ^ ^ 



r^ = F^ = U^+ J2 / dP'j,y,,.m,2l-2UP'RWR)^l-2m{P'R) ■ (57b) 
. JPX 



Lm=—oo 



In Fig. 4 we illustrated all this by the direct numerical calculation of F^ , from 
Eqs. (0) for the interaction (^), followed by a calculation of F^"^ from Eqs. (0). The RG 
solutions for F* and F* are also presented. The function F^{9) follows almost perfectly the 
Landau function (the real fixed point F*{9)), except within 29 x of 9 = 0. In the part of 
the RG trajectory px < Pr < Po (Po = 100, px = 5, T/vfKf = 0.005), not only is the ZS 
flow exponentially weak, but the central part of the ZS' flow as well (cf. Eq. (p^)). So, the 
evolution of both vertices is due mostly to the "tail" 9 > 9x of the function Y at Pr > 1. 
That is why r^{9) and F^{9) are virtually identical. Only the slowing down of the ZS' flow 
almost everywhere at Pr < 1 - except on the central part (cf. Eq. (|56|)) wherein it is always 
as strong as the other one (ZS) - results in the drastic differencesbetween the two limits 
of the four-point vertex at the fixed point. The function r^^{9) is featureless and looks like 
a corrected RPA solution. The differences between F* and F;^ {F* and F^) are negligible, 
i.e. less than 1%, only for the components n = 0,1. 

As it should be clear by now, there is no real incompatibility of the stability conditions 
with the Pauli principle, since this is a mere artefact of the ZS-ladder approximation. It is 
pointless to impose the sum rule either to in the form (^), or to F^ in the form (pSj). 
Both sums would give the value of the "uncorrelated" function r^^{9) at ^ = 0. This function 
goes smoothly from the right patch [9x, tt] towards ^ = (cf. Fig. 4) - or, equivalently, from 
the left, because of parity . Actually, it can be proved exactly, turning the arguments of 
Sec.V around, that in a stable Fermi liquid, it is impossible to obtain r^^{9 = 0) = 0, even 
by chance. Thus, there is no need for the Landau function F* to be "fine tuned" in the 
sense of the sum rule (0), since only the relation (|57| ) - between the approximate vertex 
and F^ - is an exact relationship (more precisely, asymptotically exact when T ^ 0), 
not (|3^), which relates the physical quantities F* and F*. 

In the context of our discussion at the end of Sect. |V|, notice that the cutoff Ax 
{vpAx/T ^ 1, Ax/Kp <^ 1) corresponds to the initial cutoff of the low-energy effective 
action wherein F^ is the bare interaction function (coupling) of that action. The equality of 
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the functions T-^ and illustrates the point of Sect. |T|that, at the beginning, the action's 



coupling function can be defined independently of the order in which the zero-transfer limit 
is taken. 

When the RG flow reaches the scale Ax, the contribution of the ZS' graph to the flow of 
r„ and Fn is strictly irrelevant in the RG sense, and could have been neglected in a model 
with a finite number of couplings (e.g., the theory, ID g-ology models, and so on), keeping 
only marginal terms (cf. Eqs. (^)). But, as pointed out by Shankar,EJ in the vicinity of the 
Fermi surface we are dealing with coupling functions, i.e., with an infinite set of couplings. 
Our RG solution provides a curious example of a finite deviation of the RG trajectory at the 
fixed point due to an infinite number of irrelevant terms. The right fixed point (r*(^ = 0) = 
0) canjQDt be reached if those terms are neglected, since V^^^Ox — > 0) ^ T*{9 = 0) (even at 
T = Oc^) and we would return to the problems caused by the solution F^^ (the ZS'-ladder 
approximation) discussed in Sec.V. To put it differently, neglecting those irrelevant terms 
at some part of the flow (solution ( B^)) violates the invariance of the RG trajectory at the 
point 9 = 0, expressed by Eqs. ( P3| , |26|j27] ). 

The ZS'-ladder approximation seems acceptable in the normal Fermi liquid regime with 
moderate interaction [Fn ^JS^), when the narrow-angle features of vertices revealed by the 
RG theory are not too largejHj because the forward {9 = 0) singularity has little effect on the 
first components (F* F^^, F* ^ F^ for = 0, 1 and, in the case of a weak interaction, for 
n = 2). This singularity affects mostly the higher Fourier components. So, the relationship 
(B^) is valid only for small n. It should not be used for F* [n > 2) neither directly, nor via 
the sum rule from the scattering vertex provided experimentally. For the physical vertex 



F* the sum rule 
require a large num 



is always valid, but this study indicates that its angular shape may 
Der of harmonics to adequately represent it. The existence of a finite 



solution for r^^{9) under conditions 

F^ > -1, Vn (58) 

guarantees not only finite RG solutions for F* and F*, but also the fulfillment of the ther- 
modynamic Pomeranchuk conditions (^) by F*. 

The major consequence of this study on the standard results of the Landau FLT is 
reducing the relationship (R^) between the components of the scattering vertex and the 
Landau function to the rauK of approximation and invalidating of the sum rule (pSf). The 
rest of results for normal Fermi liquids would not be affected seriously by the RG corrections. 
For example, the temperature dependence of the vertices would give a weak correction 
to the leading terms. These conclusions are neither related to the specific choice model 
considered, nor to the spatial dimension. Including spin doubles the number of vertices 
involved, changing nothing essentially. (The derivation of the RG equations with spin is 



straightforward using the A^-flavor formalism of Ref p!9|.) The differences for the case d 
are only quantitative (e.g., the type of the temperature dependence) because of different 
angular functions and solid angle integrations. 



IX. SUMMARY 



In studying the Fermi-liquid regime of interacting fermions in c? > 1 with the model of 
the ^/'^-Grassmann effective action as starting point of the analysis, one must distinguish 
between three quantities: (i) the bare interaction function of the effective action; {ii) the 
Landau interaction function; {Hi) the forward scattering vertex. We have derived the RG 
equations for the Landau channel which take into account both contributions of the ZS and 
ZS' graphs at one-loop level. The basic quantities of the Fermi liquid theory, the Landau 
function and the scattering vertex, are calculated as fixed points of the RG flow in terms of 
effective action's interaction function. 



19 



The classic derivation of Fermi liquid theory using the Bethe-Salpeter equation for the 
four-point vertex at T = is based on the approximation that the vertex irreducible in the 
direct particle-hole loop (ZS) is a regular function of its variables, neglecting the zero-angle 
singularity in the exchange loop {ZS'). This approach is equivalent to our earlier decoupled 
RG approximationllaiHy, and they are both tantamount to summation of the direct particle- 
hole ladder diagrams, wherein the Landau function stands as the bare interaction (the 
ZS*- ladder approximation) . 

One of the major deficiencies of the ZS-ladder approximation is that the antisymmetry 
of the forward scattering vertex related by the RPA-type formula to the Landau interaction 
function, is not guaranteed in the final result, and the amplitude sum rule must be imposed 
by hand on the components of the Landau function. This sum rule, not indispensable in 
the original phenomenological formulation of the Landau FLTlI, from the RG point of view 
is equivalent to fine tuning of the effective interaction. 

The strong interference of the direct and exchange processes of the particle-hole scat- 
tering near zero angle invalidates the ZS-ladder approximation in this region, resulting 
in temperature-dependent narrow-angle anomalies in the Landau function and scattering 
vertex, revealed by the RG analysis. In the present RG approach the Pauli principle is 
automatically satisfied. As follows from the RG solution, the amplitude sum rule being an 
artefact of the ZS'-ladder approximation, is not needed to respect statistics and, moreover, 
is not valid. 
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interference our results are not reliable near ^ = vr, since r(7r) = \/(0).tjll^ 



22 



FIGURES 
1 4 




FIG. 1. The three diagrams contributing to the RG flow at one- loop. 




FIG. 2. If IQ'I > 2A, the intersection (shaded) of the supports of K5 and K5 — Q' are discon- 
nected (A). If IQ'I < 2A, this intersection forms a connected area (B). Note that the RG flow is 
governed by the boundaries of this intersection, not by their interior directly. 
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FIG. 3. (A) Results of the numerical solution of the coupled RG equations. The curves labeled 
r* and F* are the forward scattering vertex and the Landau function, respectively, at tempera- 
tures T/vpKp = 0.1, 0.025, and 0.01. The narrowest central peak corresponds to the smallest 
temperature, and vice versa. (B) Approximate analytical solution of the coupled RG equations, 
for the same parameters as in (A), calculated numerically from Eq. (|49|). In both cases the initial 
cutoff was Aq = Kp. 
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FIG. 4. Comparison between the exact numerical solution of the coupled RG equations for 
T/vfKf = 0.005 (r* and F*), the intermediate values of obtained from the initial value 

U by stopping the flow at Pr = 5, and the phenomenological vertex T^^ (the result of the standard 
FLT derivations) obtained by applying the RPA solution to (F^) considered as a new initial 
point of the flow. T^^ practically coincides with F*, except in the central region. 
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